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Abstract. - We report an Inelastic Neutron Scattering (INS) study of the fully deuterated 
molecular compound Ka[Vi^ Ase042]-9D20 (V15). Due to geometrical frustration, the essential 
physics at low temperatures of the V15 cluster containing 15 coupled V 4+ (S=l/2) is determined 
by three weakly coupled spin-1/2 on a triangle. The INS spectra at low-energy allow us to 
directly determine the effective exchange coupling Jo = 0.211(2) meV within the triangle and 
the gap 2A = 0.035(2) meV between the two spin-1/2 doublets of the ground state. Results 
are discussed in terms of deviations from trigonal symmetry and Dzyaloshinskii-Moriya (DM) 
interactions. 



Introduction. - The recent discovery of quantum tunnelling of the magnetisation vector 
in magnetic cluster complexes with high-spin ground state (GS), for example Mni2-acetate 
with S = 10 has triggered a renewed interest in molecular magnets containing a 

small group of metallic ions embedded in organic or inorganic ligands. In such structures, 
each metallic cluster is magnetically isolated from its neighbours, and this allows the study of 
the collective behaviour and inherent quantum size effects within a finite number of magnetic 
atoms. Such molecular Magnets with high-spin GS prove to be extremely valuable candidates 
to study phenomena, critical at the nanoscopic scale, such as slow relaxation and quantum 
tunnelling of the magnetisation [|| , quantum coherence || and Berry phases Q . Amongst the 
existing materials, systems with low-spin GS like the polyoxovanadate K 6 [Vi5^4s6042] ' QH-zO, 
shortnamed V15 hereafter, are of the greatest interest as they are shown to display quantum 
coherence and slow relaxation of the magnetisation despite their low-spin GS and the absence 
of any magneto-crystalline energy barrier M, N] . 
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In this context, INS is an extremely valuable technique to determine the key microscopic 
parameters controlling the tunnelling mechanism in high-spin clusters such as Mni2-acetate 
H and Fes In this Letter we report the first Inelastic Neutron Scattering (INS) study 

of a fully deuterated V15 polycrystalline powder. Frustration of the exchange interactions 
plays a major role in its low-energy properties which can be accurately described by three 
weakly coupled spin-1/2. INS allows us to determine directly the effective exchange interaction 
Jo = 0.211(2) meV between the spins in the triangle and the GS splitting 2A — 0.035(2) meV. 

Magnetic properties. V15 is a polyoxovanadate molecular magnet made of 15 spins 
s=l/2 (all V 4+ ions) pi] . The symmetry is trigonal (space group R3c, a — 14.029 A a = 
79.26°). The cluster has a local C3 symmetry with the V — V distances varying between 
2.87 A and 3.05 A ; the intermolecular dipolar couplings are very small (few mK) and will 
be ignored in most situations. As first discussed in Refs. jl2| the V 4+ (S=l/2) ions are 
arranged in three layers as shown in fig. [j]a and are antiferromagnetically (AFM) coupled to 
their neighbours via oxo-bridges. Each of the two hexagons in fig. 0b contains three pairs 
of strongly coupled spins through oxo-bridge pathways between the (VO5) square pyramids. 
Three types of couplings, shown in fig. |l|b, are of importance for the magnetic properties: 
Nearest neighbour couplings J ~ 70 meV and J' ~ 13 meV and next-nearest neighbour 
couplings J" w 26 meV. In addition, each spin of the central triangle is coupled to the 
hexagons via Ji sa J' and J2 ~ J" • Crucially, the three spins S=l/2 in the triangle are not 
directly coupled to each other by oxygen bridges like in the hexagons. Two exchange pathways 
between the triangle spins should be considered: One exchange pathway occurs via the upper 
and lower hexagons with d(V — O — Vy iex — O — V) ~ 10 A while the other is mediated by 
diarsenite bridges with d(V — O — As — O — As — O — V) m 10 A . The latter coupling can 
be compared to an analogous complex, V12 |l4), where the coupling through arsenite bridges 
is J « 0.7 - 0.8 meV with d(V - O - As - O - V) « 5.3 A . The V - V distance being much 
larger in the title compound, one expects the couplings through the diarsenite bridges to be 
much smaller than 0.7 meV in V15. 

Figure. ^ shows the DC susceptibility represented as \T versus T for the fully deuterated 
version of V15. No difference in the \F could be observed with the undeuterated compound 
. The essential information from the susceptibility can be summarized as follows: At T = 
300 K the yT product does not saturate as expected for a truly paramagnetic regime thereby 
indicating that the spins are still strongly correlated. Below 100 K, the \T curve exhibits 
a plateau at around 1.1 emu.K.mol^ 1 , a value consistent with three uncoupled s=l/2 spins. 
At lower temperature, a further drop of yT suggests a S=l/2 GS. In this regime, the dimers 
in the hexagons are firmly in a singlet state (T < 0.1J) and the cluster magnetic properties 
can be accounted for by three spin-1/2 coupled by weak effective antiferromagnetic exchange 
interactions within the central triangle. An obvious starting point is then the S = 1/2 AFM 
Heisenberg model on a triangle: 

T~to = J12S1S2 + J23S2S3 + J13S1S3 : (1) 

where Si, S2, S3 denote the spin operators on sites 1,2 and 3, respectively (see fig.|l|b) and 
is the Heisenberg exchange parameter between spin i and j. We first assume the equilateral 
triangle case: Jy = Jo for all (ij) pairs. For an antiferromagnetic coupling (J > 0), the GS 
consists of two degenerate S=l/2 Kramers doublets separated from the S = 3/2 excited state 
by an energy 3Jo/2. The wave functions of the two degenerate S=l/2 Kramers doublets are 
given by *q =| 0, 5, 5, ±5 > and *f =| 1, |, i,±| > in the basis | Si 2 , S 3 , S, M s > or any 
linear combination of them. 

As a function of magnetic field the Ms = —3/2 component of the S = 3/2 spin state 
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crosses the Ms = — 1/2 level of the ground state at a critical field Hq = 3Jo/2g^B- Recently, 
Chiorescu et al. observed a jump in the low-temperature magnetisation at Hq m 2.8 T || 
which gave the first evaluation of the effective coupling: Jq — 0.211 meV. The calculated 
susceptibility using the Heisenberg model of eq. ^ and Jo = 0.211 meV matches very well 
our experimental yT versus T behaviour between 2 K and 40 K as shown in fig. |^. However 
this simple model must be refined according to several experimental results. Firstly, the 
magnetisation curve at low temperature exhibits an anomalous broadening of the jump at 
Hq Q which could not be explained by the equilateral triangular Heisenberg model and, 
more importantly, magnetisation relaxation experiments Q suggest the presence of a tunnel 
splitting 2A of the order of 9 /xeV within the S — 1/2 GS. This was ascribed to antisymmetric 
Dzyaloshinskii-Moriya (DM) terms p5| . In order to gain more insight, we have undertaken 
an INS study to directly probe the low-energy levels of this material. 

Neutron scattering experiment. We used a m *~ 4 g polycrystalline sample of fully 
deuterated V15 placed under Helium in a rectangular flat Al slab (3x5 cm) of 3 mm thick- 
ness. The INS experiments were performed on the Time-of-Flight (TOF) spectrometer NEAT 
Jl6t at the BENSC (Hahn-Meitner Institute, Germany) using cold neutrons with an incident 
wavelength of A, = 8.0 A and at the ISIS Facility (Rutherford Applcton Laboratory, United 
Kingdom) using the backscattering TOF spectrometer IRIS with a final neutron wavelength 
of Xf = 6.66 A and the pyrolytic graphite (PG002) analyzer. On NEAT, the angle between 
the incoming beam and the plane of the sample holder was 135° and the experimental energy 
resolution A_E , measured from a vanadium sample and averaged over all scattering angles, 
was 43 fieV for zero energy transfer. The energy resolution AE depends on the energy trans- 
fer hui |l7], [ll| as shown in fig. ||b. Data reduction was carried out using the program INX, 
and the fitting procedure included the energy dependence of the instrumental resolution. On 
IRIS, the experimental resolution varies from 17 fieV at the elastic position to « 28 ± 2 /ieV 
at Tvuj w ±0.3 meV. The program ICON was used for data reduction. 

The low-energy part of the scans obtained on NEAT is depicted in fig. ||a at several 
temperatures with an average momentum transfer Q av w 0.6 A _1 At T — 2 K, the 

energy loss side (left-hand side) shows a single broad peak centered at Tilo « —0.315 meV. As 
the temperature increases, a peak at the same energy appears on the energy gain side. The 
instrumental resolution at the inelastic peak positions is 32 /ieV and 55 /ieV in loss and gain, 
respectively, as shown in fig. |j^b. The observed peaks have weakly temperature-dependent 
widths with FWHM values between 57 and 76 /ieV on the loss side and between 82 and 
100 jiteV on the gain side. This is about twice the instrumental resolution, and we conclude 
that the peak contains more than one transition. In fact, the peaks have a flat-top shape in 
energy- loss suggesting that the observed broad peak is made of two unresolved INS transitions. 
Figure. ||c shows the low-energy part of the spectrum obtained on IRIS at 1.5 K and 10 K with 
with an average momentum transfer Q av w 1.1 A _1 E3|. The peaks at fko = ±0.315 meV 
correspond to those observed on NEAT but are superimposed on the thermal diffuse scattering 
of the PG002 analyzer. Again, the peaks on either side of the elastic line are broader than 
the instrumental resolution. 

The mean energy of the peak corresponds to the energy expected for the S = 1/2 «-» S = 
3/2 INS transition in the equilateral Heisenberg model {Kujq = 3Jo/2 = 0.315 meV) with 
Jo = 0.211 meV inferred from the magnetisation jump at H « 2.8 T. To account for the 
broader-than-resolution peak, we analysed the data as follows: The broad peak centered at 
0.315 meV was fitted using two Gaussians at positions ?uv = 3Jo/2 ± A. The widths are 
kept fixed at the instrumental resolution values. Only Jo, A and an overall temperature- 
dependent intensity scaling factor, It, are allowed to float. In a first approximation, It is 
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set equal for both transitions as no significant differences can be observed in their intensities 
within the experimental accuracy. The best fit is obtained for Jo = 0.2f 1 ± 0.002 meV and 
2A = 0.035 ± 0.002 meV. The results are shown as solid lines in fig. ||a. 

From these transitions alone, it is not possible to assess whether the splitting originates 
from the S=l/2 GS doublets or from the S=3/2 excited state. In both cases, the thermal 
energy, even at the lowest temperature (2 K), is much larger than 2 A « 0.4 K, and the 
change in the population factor between 2 K and higher temperatures is too small to be 
measured in the relative intensities of the two transitions fiuj^ . However, by subtracting the 
10 K data from the 1.5 K data shown in fig. ||c, it is apparent that two weak positive peaks 
occur at Tiusf = ±0.035 meV. Figure. [| shows the difference SI = Ji.sk — Jlok in the low- 
energy region | hui \ < 70 fieV. These two peaks are symmetric around the elastic line and their 
energies coincide exactly with the value of the splitting inferred from the higher energy peaks. 
Even though the accuracy of the subtracting procedure is not good enough to analyze them 
quantitatively, the fact that the transitions at fiui^ have a larger intensity at 1.5 K than at 
10 K strongly suggests that the splitting is within the two S=l/2 doublets and not within the 
S=3/2 excited state |§. 

Analysis. - To rationalise the observed splitting 2 A, the equilateral triangular Heisenberg 
model must be refined. Two types of processes can be put forward: (a) Antisymmetric 
Dzyaloshinskii-Moriya (DM) interactions [^l) , which occur as a result of spin-orbit coupling 
in non-centrosymmetric binuclear units (this is true for the edge sharing VO5 units in V15) 
and (b) deviations from purely trigonal symmetry: Removing the restriction J 12 = J13 = J23 
in eq. [l] will release part of the exchange frustration and lift the degeneracy of the S=l/2 
ground state. 

Chiorescu et al |tJ first invoked DM terms as being responsible for the splitting in the 
ground state of V15. The DM Hamiltonian may be written as 

H DM = X) 2 V0S i xS i ) ) (2) 

where Dij is the DM vector for the pair (ij) between Si and Sj along the triangle. In 
pure trigonally symmetric trinuclcar systems, Dij is identical for each bond pair and points 
parallel to the unique trigonal axis (Dij = Dz) possibly leading to a non-collinear 120°-type 
spin arrangement in the XY plane . Including DM interactions in the total Hamiltonian, 
Hi = Ho + Hdm, generates a splitting 2Ai = y/3D between the two S — 1/2 Kramers 
doublets of energy E ± — ±Ai, and the eigenfunctions are now given by specific mixture of 
*o and Hff |H, || . The observed splitting, 2 A = 0.035 meV would then lead to D w 0.02 
meV. The order of magnitude of D can be theoretically estimated from the Moriya expression 
D w (Ag/go)Jo where Ag is the deviation of the average g-factor value, g = 1.96, from 
the free-ion value go = 2.00. In our case, this would lead to a theoretical D of ~ 4.2 fJ,eV, 
a value 4.7 times smaller than the value deduced from the experimentally determined 2A. 
This discrepancy is indicative that the GS splitting is most probably not solely due to DM 
interactions. Next, we consider deviations from trigonal symmetry. Elements of non purely 
trigonal symmetry are present in the structure JllJ: Despite the C3 point symmetry of the 
complex in the average structure, disordered water molecules in the lattice can lead to local 
non-trigonal components, thereby inducing isosceles ( J12 = J and J13 = J23 = J') or scalene 
triangles ( J12 > J13 > J23 in eq. |l|) owing to the slight bond distance and angle differences. 
The GS degeneracy is then lifted with a gap 2A2 between the two Kramers doublets given by 
2A 2 = (u 2 + v 2 - uv) 1 / 2 with u = J 12 - J23 and v = J i3 - J23. No splitting of the 5 = 3/2 
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state is induced by this symmetry reduction. The lifting of the ground states degeneracy 
corresponds to a partial removal of the frustration in the equilateral triangular model. 

In the isosceles case with J > J' [|4| , the energy gap between the two Kramers doublets 
becomes 2A 2 = J — J' where is the ground state separated by 2A 2 from the tyf state. 
Using this model to analyze the INS data, we obtain J = 0.223 meV and J' = 0.188 meV, 
and the corresponding calculated yT curve between 100 K and 1.8 K is shown in fig. || 
(dashed curve). The result is not distinguishable to that of the equilateral model in this 
temperature range since T > \ J — J'\ down to T = 1.8 K. It shows, in passing, that the 
powder susceptibility alone, above T — 1.8 K (i.e. T ^> 2A), does not provide conclusive 
information on the important splitting within the 5 = 1/2 ground states p5j| . Both the DM 
interaction and a symmetry reduction of the equilateral triangle can induce a splitting in the 
ground state, and the observed 35 fj,eV gap in V15 may be a combination of both processes. 

Conclusion. - Our neutron scattering results on a deuterated powder of V15 show directly 
and unambiguously that the effective Heisenberg coupling between the 3 S=l/2 spins in the 
triangle is Jo ~ 0.211 meV and the presence of a splitting in the S=l/2 ground states of 2A « 
35/ieV. This splitting is twice as large as that reported in Ref.0] from a two-level Landau- 
Zener model analysis of the magnetisation relaxation. Our value is more accurate as it is based 
on a direct spectroscopic observation of the gap. This splitting can be attributed to deviations 
from the trigonal symmetry and/or Dzyaloshinskii-Moriya antisymmetric interactions. 
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Fig. 1 - (a) Structure of the cluster anion in [VisAsaO^] 6 ' cluster. The 15 V 4+ (dark solid spheres) 
form two outer distorted hexagonal layers and an inner triangular layer, (b) Schematic represen- 
tation of the AFM exchange couplings between V 4+ ions according to the type of oxygen bridges 
involved in the exchange paths (distance and angle). The magnetic Hamiltonian is defined as 
Ti — J^<j > JijSiSj with < i,j > pairs as indicated. Si, S2 and S3 in the triangle are loosely 
connected through the hexagons and diarsenite bridges. 
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Isosceles Heisenberg 
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Fig. 2 - DC SQUID magnetic susceptibility obtained on a 160 mg powder sample of deuterated V15 
at an applied external field H — 0.1 T. The powder averaged Lande g-factor is set to 1.96 from EPR 
data: g a = Qb = 1.95, g c = 1.98 jl^, |l3). Solid and dashed lines are calculated curves from eq.[l] for 
an equilateral triangle (J12 = J23 = J13 = 0.211 meV) and an isosceles triangle (J12 = J23 = 0.223 
meV and J12 = 0.188 meV), respectively (see text). 
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Fig. 3 - (a) INS spectra of fully deuterated Vis obtained on NEAT (HMI) at A* = 8.0 A with detector 
angles 23° < 29 < 11° . Data were collected at several temperatures and corrected for the background 
(empty sample holder) and detector efficiency (Vanadium metal reference). The solid lines are fits to 
the data as explained in the text, (b) Instrumental resolution AE as a function of the energy transfer 
Tiu> calculated for NEAT. The calculation of the energy dependence of A_E is performed according to 
NEAT specifications |l7]] and rescaled by hand to match the measured elastic energy resolution A£o- 
(c) INS spectra at 1.5 K and 10 K obtained on IRIS (ISIS) with Xf = 6.66 A _1 and for detector 
angles 28.5° < 29 < 128°. Data were corrected similarly as for NEAT. The T = 1.5 K Data after 
subtraction of the thermal diffuse scattering is shown a solid lines. 
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Fig. 4 - Difference between the 1.5 K and the 10 K IRIS data. The arrows indicate the positive peaks 
discussed in the text. 



